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ON  SOME  JOINT  PROBABILITIES  USEFUL  IN  MIXED  ACCEPTANCE  SAMPLING 

by 

E.  G.  Schilling1 
and 

H.  F.  Dodge 
1.  INTRODUCTION 


Mixed  acceptance  sampling  plans,  i.e.  sampling  schemes  combining  constituent 
variables  and  attributes  results  in  a  double  sampling  procedure,  provide  con¬ 
siderable  advantage  over  more  conventional  plans  in  certain  applications.  This 
has  been  pointed  out  by  Dodge  Tlj  and  Bowker  and  Goode  [2],  among  others.  Such 
plans  have  not,  however,  found  application  in  proportion  to  their  merit,  partly 
because  methods  of  computing  the  probabilities  necessary  for  the  evaluation  of 
the  operating  characteristics  of  such  plans  far  all  attributes  acceptance  nuirbers 
have  not  been  developed. 

Of  particular  importance  is  the  case  involving  a  single  specification  limit, 
known  standard  deviation,  assuming  an  underlying  normal  distribution  of  product, 
when  the  results  of  the  variables  and  attributes  constituents  of  the  plan  are 
not  kept  independent.  At  present  the  only  detailed,  investigation  in  this  area 
appears  to  be  the  pioneering  work  of  Gregory  and  Resnikoff  [3j*  Savage  TU]  points 
out  that  Gregory  and  Resnikoff  were  able  to  find  feasible  computing  techniques  for  the 
necessary  probabilities  for  this  case  only  when  c  ■  0,  a  considerable  restriction. 

This  report  includes  an  alternative  derivation  of  the  necessary  probabilities  when 
c  ■  0  and  then  presents  a  derivation  of  the  probabilities  for  c  >  0,  together  with 

^The  material  in  this  report  is  based  in  part  on  work  done  in  preparation  of 
a  doctoral  dissertation  at  Rutgers  -  The  State  University. 

Distribution  of  this  document  is  unlimited. 
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a  scheme  for  readily  computing  these  probabilities. 

1.2.  The  Mixed  Acceptance  Sampling  Plan 

Mixed  acceptance  sampling  plans  are  of  two  types,  •‘independent*  and  "dependent 
Independent  plans  are  those  in  which  the  results  of  the  variables  and  attributes 
components  are  kept  independent.  Such  plans  have  been  discussed  by  Bowker  and 
Ooode  [1],  Oregory  and  Resnikoff  [3],  and  Schilling  [*>]•  The  values  of  probability 
of  acceptance,  Pft,  of  such  plans  are  readily  determined  through  the  use  of  existing 
tables  and  procedures.  Probably  of  more  importance,  both  theoretically  and  practi¬ 
cally,  are  dependent  plans,  those  in  which  the  variables  and  attributes  results  of 
the  procedure  are  combined  in  such  a  way  as  to  make  them  dependent.  Attention 
will  be  directed  here  to  dependent  plans  involving  a  single  specification  limit, 
known  standard  deviation,  and  the  assumption  of  an  underlying  normal  distribution 
of  product.  For  reasons  of  symmetry,  only  an  upper  specification  limit  will  be 
considered. 

Let: 

n^  -  the  first  sample  size. 

A  ■  variables  acceptance  limit  such  that  if  the  sample  mean, 

-  2 

x  £  A  the  lot  is  accepted  . 

n2  *  second  sample  size, 
c  *  attributes  acceptance  number. 

There  are  a  number  of  procedures  for  carrying  out  a  dependent  mixed  plan.  A 
typical  one  is  that  originally  proposed  by  Bowker  and  Goode  Cl]  which  can  be 
summarized  as  follows: 

2 

Of  the  several  methods  of  specifying  the  variables  constituent  of  known-standard- 
deviation  (a)  variables  plans,  designation  by  sample  size  (n^)  and  acceptance 

limit  on  the  sample  average  (A),  simplifies  the  notation  somewhat.  Note  that 
A  *  U-ko  for  upper  specification  limit  U  and  customary  variables  acceptance 
factor  k$  also  that  A  ■  (K  -  k)  in  the  notation  of  Gregory  and  Resnikoff  [31# 
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1.  Determine  parameters  of  plan:  n^,  A,  r^,  c. 

2.  Take  a  random  sample  of  size  from  the  lot* 

3.  Determine  the  sample  average  (x)  of  this  first  sample  and, 

a.  If  5c  <  A  accept  the  lot  and  cease  inspection. 

b.  If  5c  >  A  continue  inspection  as  in  U* 

U,  Take  a  second  sample  of  size  from  the  lot  and  determine 

the  number  of  defectives  m  contained  in  the  total  combined 
sample  of  and, 

a*  If  m  <  c  accept  the  lot  and  cease  inspection* 
b.  If  m  >  c  reject  the  lot  and  cease  inspection. 

Clearly,  since  the  combined  sample  is  subjected  to  the  attributes  test,  the  vari¬ 
ables  and  attributes  aspects  of  the  plan  are  not  independent.  The  probability 
of  acceptance  under  this  procedure  becomes: 

c  c— 1 

P.  -  P(x  *  A)  ♦  r  T,  P^ (i,x  >  A)  PU;nJ 
a  1-0  J-0  "l 


where 

P(5c  s  l)  •  probability  that  the  sample  mean  does  not  exceed  the 
acceptance  limit  A 

pCjjng)  ■  probability  of  J  defectives  in  a  sample  of  size  rx^ 

P  (i,5c  >  A)  ■  Joint  probability  of  5c  >  A  and  i  defectives  in 
i  a  sample  of  size  n. * 


Examination  of  the  formula  for  the  probability  of  acceptance,  P^,  will  in¬ 
dicate  that  all  probabilities  contained  therein  have  been  tabulated  or  can  be 

calculated  by  elementary  methods,  with  the  exception  of  P  (i,5c  >  A)*  The  same 

°1 

would  be  true  for  modified  forms  of  the  above  procedure,  say  one  having  two  ac¬ 
ceptance  numbers,  and  Cg,  for  tho  2  samples.  Thus,  this  joint  probability  is 
a  very  Important  factor  in  determining  the  properties  of  mixed  sampling  plans 
generally.  It  provides  the  key,  and  can  be  used  directly,  for  determining  the 
operating  characteristic  curves  of  mixed  plans  of  the  general  type  to  be  considered 
here  —  for  single  specif ication  limit,  and  known  standard  deviation. 
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2.  METHOD  OF  EVAHJATIl©  P„(i,x  >  k) 

n 

2.1,  Approach 

This  section  will  present  a  method  of  evaluating  Fn(i,x  >  A),  the  probability 
that  a  sample  of  size  n  will  have  exactly  i  defectives  and  a  sample  mean  (x) 
greater  than  some  limit  A,  First  to  be  considered  will  be  the  case  in  which  i  ■  0, 
This  will  be  followed  by  the  case  in  which  1  >  0,  The  method  of  evaluation  used 
was  suggested  from  the  results  of  Gregory  and  Resnikoff  [3]  —  the  solution  fcr 

3 

this  probability  when  i  ■  0  is  the  same  a3  that  of  their  exposition  , 

In  determining  these  probabilities,  we  will  let  pn(i,x)  be  the  joint  prob¬ 
ability  density  of  i  and  x  in  a  sample  of  n.  Then,  given  a  particular  value 
i  •  ig,  integration  of  this  density  over  an  appropriate  region  for  x  will  give 
Pn(i,x  >  A)  if  the  function  is  integrable  and  the  integral  exists. 

2.2,  Evaluation  of  Pn(0,x  >  A) 

Without  loss  of  generality,  since  the  z  transformation  will  always  be 
appropriate,  consider  random  samples  of  n  from  the  standard  normal  distribution* 
5(0,1),  Choose  z^  so  that: 


e  dt  ■  p 


Thus,  8^  is  an  upper  limit  such  that  only  lOOp  percent  of  the  population  exceeds 
the  said  limit.  The  probability  density  for  samples  with  mean  z  and  no  obser¬ 
vations  greater  than  z^  is  desired,  that  is,  samples  with  mean  z  such  that  if  the 
observations  are  ordered! 


3 

The  Gregory  and  Resnikoff  evaluation  of  a  similar  probability  for  i  >  0  appears 
to  be  subject  to  question  due  to  the  particular  restrictions  placed  on  the 
region  of  integration. 


TECHNICAL  REPORT  ND,  N-26 
Issue  1  -  Page  5 


*(1)  *  *(2)  <  *  •  *  <  *(n)  <  * 

Now  the  probability  that  none  of  the  observations  in  a  random  sample  from 
♦(0,1)  is  a  distance  further  from  the  sample  mean  than  the  limit  z^  can  be  found 
from  the  distribution  of  the  deviation  of  z ^  from  the  sample  mean,  uR  » 
reouiring  that  (z^-z)  <  (Zy-z).  This  distribution  has  been  tabulated  both  by 
Nair  [6'J  and  Grubbs  [7]  as  Fn(u)  for  *(0,1)  where  Fn(u)  ■  P(un  <  u)  for  a  sample 
of  size  n. 

It  can  be  shown  that  the  extreme  deviate  from  the  sample  mean,  un,  and  the 
sample  mean  »  are  Independent.  Consider  the  following  theorem! 

Th.  Let  X^Xg,...,^  denote  a  random  sample  from  a 

distribution  having  a  p.d.f,  f(xj0)  y  <  P  <  6.  Let 
-  u^(X^,l2>*..#Xn)  be  a  sufficient  statistic  for 

0  and  let  the  p.d.f.  g^Cy^jQ^Y  <0<6,  ofT^be  complete. 

Let  Z  ■  uClpXg,...,!,^)  be  any  other  statistic  (not  a 

function  of  alone).  If  the  distribution  of  Z  does  not 

depend  upon  0,  then  Z  is  stochastically  independent  of  the 

sufficient  statistic 

-  —  2 
Here  *  a  irfiile  Z  -  un.  Also  z  is,  for  every  a  >  0,  a  sufficient  statistic 

for  \i9 <  ja  <  ®,  and  the  density  function  of  the  sufficient  statistic  *  is 

complete.  The  cumulative  distribution  function  of  u  is 
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which  does  not  contain  and  hence  uR  ■  (z^-i)  is  stochastically  independent 
of  the  sufficient  statistic  z. 

So  when  i  ■  0,  the  joint  density  for  z  in  the  range  -®  <  z  <  is 

pn(0,¥)  -  p(z)  P(*(n)“*  <  zn-*)  "  P(*)  Fn(zu-I). 

But  Pn(0,¥)  is  Riemann  integrable  in  the  region  -®  <  z  *  since  it  is  the  product 
of  two  continuous  functions  in  the  interval*  For  some  value  z^  such  that 

-'*A<  V 


P»(0'*  >  *A> 


J 


rS2 

*0  "  V 

/  ^  .  F„(v*> 


n  U 


and  the  right  hand  side  can  be  evaluated  by  numerical  integration  procedures. 

2.3.  Evaluation  of  Pn(i,x  >  A),i  >  0 

Again  without  loss  of  generality,  consider  random  samples  from  the  standard 
normal  distribution  $(0,l),  The  probability  density  for  samples  with  mean  z  and 
i  observations  greater  than  z^  is  desired,  that  is,  samples  with  mean  z  such  that 
if  the  observations  are  ordered: 

-  <  ‘(1)  ‘  l(2)  4  —  *  *(n-i)  4  *U  4  z(n-i+l)  4  ”*  4  *<n)  4  ”  * 

Suppose  the  sample  to  be  selected  sequentially.  Further,  suppose  that  the  sequence 
of  observations  givir.g  rise  to  such  a  sample  was  in  the  order  indicated,  that  is, 
the  first  (n-i)  observations  less  than  or  equal  to  z^  and  the  last  (i)  observations 
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greater  than  s^,  Then  the  sample  could  be  divided  Into  two  subsamples,  the  first 
(n-i)  values  <  and  the  second  (i)  values  >  If  these  subsanples  are  regarded 
separately,  the  means  of  the  two  sang>les  would  be  related  to  the  overall  nee*  of 
the  grand  sample  ast 

(n-i)i^  +  (i)z2  -  (n)¥ 


while  clearly 


-  <  *1 4  *0 

*0  <  *2  <  -• 

Now  to  have  a  grand  mean  ¥  and  1  defectives,  there  must  be  associated  with 
?  ,  no  sample  value  >  By  and  with  *2  no  sample  value  <  By.  But  since  the  grand 
sample  has  been  divided  into  two  subsamples,  each  of  them  can  be  treated  in 
exactly  the  manner  of  the  above  section  for  Pn(0,*).  So  for 

Vi(°«Ti>  *  WvV' 

And  treating  the  lower  limit  on  the  second  sample  as  If  It  were  an  upper  limit! 


*  P(*2^i^*2"*U^ 

since  Fn(u)  ■  P(un  <  u)  •  P(-u  s  u^)* 

Note  that  in  the  above  section  it  was  shown  that  i*  independent  of  a 

for  the  son  pie.  Since  each  of  the  subsamples  is  being  treated  as  an  entity  in 

Itself,  the  magnitude  of  7,  will  not  affect  F  (u).  This,  of  course,  must  be  the 

x  n 

case  for  the  statistic  (■  ■  '  )  to  be  useful  in  testing  outliers,  its  principal 

c 

application. 
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For  a  grand  mean  *  in  the  range  <  z  <  ®  and  z^  and  >2  related  as  shown 
above,  the  joint  density  is  such  that 

•pn(i,z)  -  (J)  Pn.1(0»*1)p1(i#*2^ 

"  (i)P^i)Fn.j.(*u“*i)P^2^Fi^*2“aU^ 

since  there  are  (“)  possible  combinations  of  the  observations  to  give  n-i  <  U  and 

i  >  U. 

Now  pn(i,z)  certainly  does  not  define  a  continuous  surface.  However,  for 
a  fixed  value  of  i  ■  1^,  it  becomes  a  product  of  continuous  functions  and  is 
therefore  Rlemann  integrable  over  an  appropriate  finite  range.  Furthermore,  it  is 
dominated  by  an  appropriate  multiple  of  the  normal  density  and  therefore  the  im¬ 
proper  integral  exists.  Again  for  any  specified  value  z^,  <  \  <  ",  Pn(i#») 

must  be  integrated  over  a  suitaole  region  to  obtain  the  desired  probability, 

Pn(i,z  >  a^).  Given  z,  if  z 2  is  allowed  to  vary  so  that  zy  <  a2  <  »,  then 
will  be  determined  by  the  relationship* 


I  -  xs. 

ST”  <  ?1  ‘  V 


For  Pn(i,z  >  however,  z  is  never  <  »^,  z^  being  the  greatest  lower  bound 
for  the  interval, 
so 


na.-iZp 

“ i  ”  <?1S  ®A* 


Hence,  for  i  »  1^  >  0  and  appropriate  regions  of  integration  for  z^  and  *2, 
and  PSg  respectively* 


l 
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V1'*  >  zk]  ’  J  I  ( pf*2 >7i ( *2-“u)p C*l)Frwi< ,U-5l)d*l<1*2 


FG?2  Rz^ 


/S3 


(?) J  l  %•  Fi<w • 

Rs2  R«1 


-  ^  “7  [l*2Z+<'*-i)*12] 

(?)  I  J  _  ^  •  VW^vW* 

*U  nZA“iz2 

"~H=T“ 


Which  can  be  integrated  numerically  as  will  be  shown* 

2.U*  Summary 

The  probability  Pn(i,T  >  a^)  can  be  evaluated  by  use  of  the  formula  ex 


n* 

*U  "  T 

Pn(0,5>.A)-J 


Fn(Vi)'ff 


i  -  0 


and,  for  i  ■  iQ  >  0 


Pn(1«*  >  **) 


■o 


(?)/  J 

*u  mCsz 


/i(n-i)  . 

~as —  • 


-  |  [l^2*(n-i)*12  ] 


h^-'J  Fn-i(V1‘l)  S1S2  •  1  >  0 

Pn(0,z  >  z^)  and  Pn(i,a  >  z^),  then,  represent  values  of  the  probability  desired 
for  samples  taken  from  a  standard  normal  universe.  Using  the  s  transformation 


for  a  particular  i  ■  i^x 
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Pn(i,5  >  A)  -  Pn(i,s  -  >  sA  - 

and  conversely 

Pn(i,s  >  »A)  -  Pn(i,x  -  ab  ♦  n  >  A  -  *Ao  ♦  u) 

3.  COMHJTATION  OF  Pn(i,x  >  A) 

3.1.  Scheme  for  Computation 
The  computation  of 


P„<°,5 


V 


*o 

f  / n 

J  75T 


-m2 


Fn(VI)<!l* 


can  be  accomplished  by  the  application  of  appropriate  methods  of  numerical  Inte¬ 
gration  to  the  integrand  between  the  limits  indicated.  Values  of  F^a^-a)  may 

be  obtained  by  interpolation  from  tha  tables  of  F  (u)  of  both  Nair  T61  and 

n 

Orubbs  [71  to  arrive  at  the  functional  values. 

For  a  particular  i  -  i0,  the  computation  of 


"  *D  _  -  j  [llj^tn-l)^2  ] 

pn(u  >  v  *  <;>  I  J 

»o  Mr1?2 

n-i 


is  greatly  simplified  by  rewriting  the  integral  as 
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Pn(i,S  >  *A) 


112 
"  7  *2 


(i}  1 75T  *  fi(?2"*u) 

*o 


*u 

/r£i  f 

7sr  J 


_(n-i)*12 

*  WVW 


ml'U2 

n-1 


and  noting  that 


»  -li2 

>  -sa-1) 


-(n-i)I  2 


J  7p‘  ‘  WV*>* 


n*  -is_ 

A  2 

n-i’~ ' 


Whan  calculated  for  a  sample  of  slxe  (n-i).  This  Is  exactly  the  expression 
which  appears  In  brackets  In  the  formula  for  P(l,e  >  «^)  a  bore. 

So, 

.  1  a  2 

“7  *2 

%  •  yi<VV  [Pn-i(0'Il  > 

*0 


P(i,7  >  «,)  -  (J)  { 


Treating  the  Integrand  as  one  function  for  purpose  of  numerical  integration  it  is 

possible  to  interpolate  for  values  of  P^CO,*^  >  — gjy . )  once  tables  of  this 

probability  are  available.  In  the  same  way  ?A(*2“*u)  **  obtained  from  existing 

tables,  where  available,  to  arrive  at  values  of  the  integrand. 

In  evaluating  the  integral  expression  for  Pn(i,*  >  *A)  it  is  worthwhile  to 
note  that  the  values  of  i  for  which  F^Ca^s^)  is  calculated  will  generally  be  sdmII. 
Bat  F^(u)  “  1  for  all  values  of  u,  since,  for  n  ■  1,  it  represents  P(x^-x^  <  u)  for 
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u  *  0.  Further,  for  i  -  2* 


F2(u) 


-  2 


u 

n  f 


-  i(2*2) 


7S 


dx 


which  is  simply  the  error  function  for  which  standard  approximations  are  available. 
These  relationships  simplify  the  computation  of  Pn(i,z  >  z^)  still  further. 

5.  CONCLUSION 


This  report  presents  a  method  for  evaluating  joint  Drobabilitles  necessary 
to  determine  the  operating  characteristics  of  dependent  mixed  plans  for  the  case 
of  single  specification  limit,  known  standard  deviation,  assuming  an  underlying 
normal  distribution.  It  also  provides  a  technique  for  facilitating  the  computation 
of  these  joint  probabilities.  3uch  values  will  be  of  substantial  aid  to  any  com¬ 
prehensive  studies  that  may  be  made  to  determine  the  best  composite  structure  and 
the  relative  merits  of  denendent  mixed  plans,  or  to  explore  the  possibility  of 
developing  suitable  sets  of  mixed  plans  for  practical  application.  Tables  of 
values  for  a  number  of  small  sample  sizes  are  being  developed. 
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